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Abstract 



Bannonl [20021 . J. Atmos. Sci. 59: 1967-1982] developed equations governing the 



motion of moist air. These equations include a term for a "reactive motion" that arises 
when water vapor condenses and droplets begin to fall; according to this term the 
remaining gas moves upwards so as to conserve momentum. Here we show that this 
reactive motion is based on a misunderstanding of the conservation of momentum in the 
presence of a gravitational field. Other suggested means to derive the same expression 
are shown to be based on an incorrect application of Newton's second law. We conclude 
that the assumed reactive motion is not compatible with the fundamental physical 
principles: it does not exist. 



1 Introduction 

The foundations of meteorology and climate sciences lie in the physical prin- 
ciples that govern atmospheric behavior. Basic physical principles, such as the 
conservation of momentum, must be correctly understood and applied. Here 
we identify and explain a key error in one published fundamental equation 
go verning atmosp heric motion that has persisted in the literature for a decade. 

investigated how the Navier-Stokes equations that determine 
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velocity u a of moist air 

Pm (-t^ + (u • V)u a J = -Vp + V • a + p m g 

(here p m = p a -\- p v , p a and p v are moist air, dry air and water vapor density 
(kg m~ 3 ), p is air pressure, o is the viscous stress tensor, g is acceleration of 
gravity) are modified by the presence of condensate particles in the atmosphere. 
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Bannon 2002 concluded that condensation causes an additional term F„ to 



appear on the right-hand side of Eq. 



F a - pig + vipv 



(2) 



Here p\ = Nffn is the density of condensate particles in the air, Ni is the 
concentration of condensate particles in the air, fn is their mean mass, p v < 
(kg m~ 3 s" 1 ) is condensation rate and v/ is the mean droplet velocity relative 
to the air. We shall consider only liquid droplets explicitly without losing the 
generality of our conclusions. 

The secon d term in Eq. (El) describes an upward acceleration of air upon 
condensation. 



Bannon 



20021 . p. 1972] reasoned the need for this term 



active motion": with the droplets acquiring a downward velocity, the remaining 
air must gain an upward velocity so that total momentum (air plus droplets) is 
conserved. Here we shall show that in the presence of a gravitational field this 
supposed reactive motion contradicts physical principles and does not exist. 



2 Misunderstanding the conservation of momentum 



Consider the argument for the reactive motion term in Eq. ([2]). iBannonl [20021 . 
p. 1972] discussed a body of mass M that is subjected to no external forces. 
For an example of such an idealized system, we can take a rocket that moves 
with velocity u a in the absence of gravity. The rocket's motion is reactive: 
it expels some of its parts that move at a different velocity u;. The mass of 
the rocket diminishes at a rate dM/dt < 0. Neglecting the acceleration of the 
expelled fuel/exhaust, dui/dt = (cf. the forthcoming Eq. (JED), we can write 
the momentum conservation equation for the rocket and all of its parts as 

d(Mu a ) dM 

M^ = v^, (4) 
dt dt ' w 

where = u; — u a . Thus, acceleration of the rocket is proportional to the rate 
at which it loses mass and to the difference v/ between the rocket velocity and 
the velocity of its expelled parts. Equation (Jlj) describes the reactive motion. Its 
right-hand part is fully analogous to the second term i n Eq. (j2 ) whil e Eq. 



2002 



is equivalent to Eq. (B5.16) ("B" refers to equations of iBannon 

However, unlike our hypothetical rocket, atmospheric air does not exist in 
a space devoid of external force - it is subjected to gravity. The conservation 
of momentum applicable to our rocket +exhaust system is not applicable to 
the air alone. The Earth's gravitational field is not a detail here: it is central to 
the process. Unlike the fuel/exhaust that is expelled from the rocket, the newly 
formed droplets accelerate freely in a uniform downward direction not under 
some internal process of physical chemistry (e.g., fuel combustion), but under 
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the pervasive influence of the Earth's gravity. Certainly the momentum of the 
system is conserved but here "the system" includes the entire bulk of the planet 
Earth. Based on these considerations we conclude that the proposed reactive 
motion of moist air based on the supposed conservation of momentum of the 
isolated system (air+droplets) has no physical basis. It does not exist. 



3 Examining the derivations of iBannonl |2002|| 



How did the "reactive motion" term arise ? IBannonl [20021 ] . following the general 
approach outlined by ICrowe et al.l |1998l | , based his derivation on the following 
key equation (B2.2): 



— + V 

dt 



(x%) = x, 



(5) 



This equation represents a general conservation relationship for an arbi- 
tra ry property y that is transported with velocity u x . Here x was described 
by IBannonl [20021 ] as the rate at which property x is being created within the 
considered unit volume. A problem with Eq. (jSj) is that x remains undefined. 
For any property x-> one must refer to some relationship (here a physical law) 
to determine x- 



We note that Eq. © or (B2.2) is equivalent t o Eq. (3.1 . 16) o f 



1967 1 with x^^P an d X = Q m the notations of iBatchelor 



1967 



Batchelor 



Batchelor 



19671 . p. 136] specifically mentioned the problem that x depends on the nature 



of the considered quantity x such that no general relationship for it exists. 

Intending to apply Eq for the to tal momentum of a unit atmospheric 
volume, x = Pm^-a + Pi**-!-, IBannonl |2002| did not specify, discuss or justify any 
such relationship. He disregarded the problem mentioned by others [Batchelor! . 
1967fl. Instead, he postulated that Eq. (J5J) written for x = Pm u a + Pi^i is equiv- 
alent to the following relationship, see (B5.2) and (B5.7): 



Du a Dui 

Pm—RT + Pl-RT = ~VP + V • O + p m g + pig + Vip v . 



Dt 



Dti 



D d . 
D-t = di + {Ua 



Here u; is the mean velocity of droplets in a unit atmospheric volume. O 
;hat d roplets fall at a nearly constant velocity for most of their lifetime 



(6) 
(7) 

Dserving 



Bannon 



20021 . (B5.11a)] proposed that 



Put 

Dt, 



0. 



Using Eq. (j8J) in Eq. (J6J) and comparing Eqs. (jbj) and (JTJ) we obtain: 

DUr. 



P, 



Dt 



- = -Vp + V • a + p m g + F , F a - pig + vip\ 



(9) 
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Equation (Ej) is equivalent to Eq. (B5.18). (In the notations of (B5.18) we have 
Pv = -A^cond and F a = p a u = pig - p a u ph ^c = pig + vip v , see (B5.18) and 
(B 5.7b).) 



Bannonl |2002| never discussed the problem of x being undefined but rather 



based his derivations on an incorrect postulate. 



4 Res olving t 



and Bannon 



l e con tradiction between lOoyamal |200ll | 
|20Q2l | 



Our evaluations help illuminat e a relat e d top ic: t he contradiction between the 
expressions for F „ obtained bv IBannonl (20021 ] and lOoyamal |2001| . The expres- 



sion obtained by lOoyama 



2001 



does not contain the reactive motion term. 

the 



While the contradiction is recognised [Bannonl . 120021 . ICotton et al. 
implications remain uncertain. 



2011 



Ooyamal |2001| assumed that in the presence of condensate the equation for 
the motion of moist air contains an additional term F a that describes the rate 
of momentum transfer from the droplets to the air: 

Du n 



P, 



Dt 



-Vp + V -a + p m g + F c 



(10) 



ikewise, the equation for the motion of droplets contains, following lOoyama 



2001 , an additional term F; that describes the rate of momentum transfer 



from the air to the droplets^: 



(u) 



Without addressing the nature of this momentum exchange, lOoyamal [2001 
assumed that they represent internal interactions within the considered atmo 
spheric volume and cancel out when Eqs. ( TTT1 ) and ( TTUi ) are summeco: 



= -x i- 



Summing Eqs. ( TTT1 ) and ( TTUI ) in the view of ( TT21 yields 



P, 



Du a Dui _ _ 

'"DT + Pl 15ti = + V • CT + p m g + pig. 



(12) 



(13) 



As one can see, Eq. ( TTcfl ) differs from Eq. of IBannonl |2002| in that in the 
for mer the reactiv e motion term is absent. If lOoyama |200l| had assumed, as 
did IBannonl |2002| , that Dui/Dti = (ED, his expression for F a would have 
been F a = pig 7 cf. ([JU]) and (H 



1 See Section ^ for a discussion of this equation. 

2 We plan to examine this equ ation elsewhere - its validity does not impact our conclusions here. 

3 To be accurate we note that Oovamal [200 1 1 made a slightly different assumption, Dvi/Dti = 0, which 
was presumed to reflect no change of terminal velocity along the droplet path. Accordingly, his expression 
for F a in (TTU)) became F a = pig — piDu a / Dti. This detail is not however relevant to the reactive motion 
problem. 
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Th us w e can se e that the controversy between the derivations of lOoyama 
200 1| and lBannonl |2002| is resolved. It results from the incorrect equation (ED 
of Bannonl 20021. 



5 An alternative derivation of the reactive motion term 

One anonymous referee^ suggested a procedure to derive Eq. (J2D without the 
problematic Eq. ([5]). Here we will show that this proposed derivation is also 
based on a physical error and cannot be justified. The referee started by noting 
that Newton's second law for a particle of mass M and velocity u has the form: 

where f ex t represents all external forces acting on the particle. Observing that in 
a continuous medium the momentum is defined as J pudr, where p is density 
and dr is an elementary volume of moving fluid, the referee presumed that 
Newton's second law, by analogy with (TH1). can be written as: 

— J pudr = f ex t = J~F ext dr, (15) 

where F ext is the external force per unit atmospheric air volume. 
Noting that 



UpudT= /pgW /"u(g + p(V.u))rfr, (16) 



see Eq. (3.1.11) of lBatchelorl [1967| . taking into account the continuity equation 

d Dp 

j t {pdr) = pdr or — + p(V • u) = p (17) 

and assuming that Eqs. ( TT5T ) and ( TMj ) are valid for any dr, they derived the 
following equation from Eqs. ( 1TBT) and (TT51): 

Pj^T = F ex t ~ pu. (18) 

Writing Eq. ( TTH1 ) for the air (u = u a , p = p m ) and droplets (u = u/, p = pi) 
yielded: 

pi—^ = Pig + F fl + p v ui, (20) 

Here Fj a and are the friction forces exerted by the droplets on the air 
and vice versa, respectively. These forces exist independent of condensation. 



'See details here. 
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By summing Eqs. ( TT§1 ) and ( 1201 ). taking into account tha t v; = u/ — u „ and 



Bannon 



2002 



assuming that F/ a = — F/j the referee obtained Eq. (ED of 

The error here is that Newton's second law in the form (Il4l ) only applies for 
constant mass systems where dM/dt = 0. In the general case of bodies with 
variable mass the correct form of Newton's second law is 

Similarly, in the dynamic Euler and Navier-Stokes equations (OQ) density p in 
the pDu/Dt term is not under the sign of differentiation irrespective of whether 
p is variable (including possible changes because of condensation) or constant. 
The fundamental dynamic equations based on Newton's second law are the 
equations on velocity not momentum. For momentum, there are no dynamic 
equations, that is, no equations defining \ m Eq. (ED for x representing mo- 
mentum. For momentum, there is only the law of momentum conservation. 

In the view of Eq. ( |2~TT) . equation ( TT51 ) is only valid when the mass of the 
fluid with velocity u is conserved, i.e. when pdr is constant over time and p in 
Eq. ( ITTjl is equal to zero. In this case Eq. ( TTH1 ) corresponds to the fundamental 
Euler or Navier-Stokes equation. When p ^ 0, Eq. (JT5T) as well as the subsequent 
derivations are invalid. 

The incorrect application of Newton's second law by using the form (TT4l in 
cases where mass is not fixed is surprisingly common. It has been discussed 
in journal a rticles as well as in textbo oks: for details we refer the reader, for 
example, to Plastino and Muzziol 1992 ] and references therein. As a brief illus- 
tration of the absurd unphysical results that follow from Eq. ( TT4l ) when mass 
is changing, consider a block of ice of mass M that is melting in a warm room, 
dM/dt < 0, in the absence of net external forces, f ext = 0: 

, ^du dM 

M— = -u— . 22 

dt dt v ; 

Since velocity u = 0, we have from (J22D Mdu/dt = 0: no acceleration, the 
ice block remains motionless on the table. Consider now the same piece of 
ice travelling in a train that moves with constant velocity u ^ 0. Then we 
obtain from (TT4l that Mdu/dt = —udM/dt ^ 0. This means that the ice 
block accelerates within the train because of melting. This absurd result (non- 
invariance of acceleration in different inertial frames of reference) illustrates the 
flaw ed logic behind applying Eq. f lT4l to systems with changing mass. We note 
that iBannonl |2002| followed the same incorrect logic. His unnumbered equation 



in the right column on p. 1972 is identical to Eq. ( |2"2~p . He also confused Eq. 
which has the same velocity u on both sides of the equation, with Eq. (jjj) for 
reactive motion, where the velocities are different. 

Another example of unphysical results stemming from Eq. ( TT5| ) at p ^ 
is given by Eq. (l20l) . Indeed, this equation, derived by the referee, contradicts 
the commonly used equation of motion for droplets, where the p v \ii term is 
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normally abse nt. S pecifically. Eg. (12Uj) contradi cts Eq. (B5.8a) borrowed by 
Bannonl |2002l | from IPruppacher and Klettl [1997 . 



6 The motion of condensate particles 

We could stop our analysis here. However, it may be instruc tive to analyze what 



considerations could motivate the derivations like that of iBannonl |2002l | that 
start from formulating the conservation equations (J5J) for total momentum. Such 
derivations could be undertaken to determine the links between the movement 
of air and of droplets and to describe the exchange of momentum between them. 
Indeed, the continuity equations for water vapor and liquid, 

^ + V • (p v u a ) = p v , (23) 

^ + V • (pm) = -p v: (24) 

seemingly tell us that condensation turns water vapor that moves with velocity 
u a into liquid droplets that move with mean velocity u; ^ u a . As a consequence, 
the change of total momentum appears to depend on condensation rate p v and 
on the difference in velocities v/ = u; — u a . However, as we have discussed, such 
an exchange of momentum does not in fact take place. 

Motion of any particular droplet is described by the linear equations of 
Newton's second law; it is independent of the presence of other droplets. For 
this reason there cannot be any non-linear dynamic equation for the mean 
velocity of droplets in a unit volume. (In contrast, for air such a non- linear 
equation does exist: it is the Euler equation.) With no dynamic constraints 
on u;, the continuity equation for droplets ( 1241 that includes the undefined 
u; becomes meaningless. Any physical inferences drawn from this equation are 
potentially misleading. We will now prove that the non-linear dynamic equation 
( TTT1 ) for condensate particles is incorrect. 

Consider a flux of macroscopic particles (raindrops, hailstones, etc.) that do 
not interact with each other contained within a fixed volume. We can define 
the mean velocity u; and mean acceleration a; of particles as 

u, = £=i!* a, = £=1* a„ = ^. (25) 
n n at 

Here n is the number of particles in an arbitrary small volume V enclosing the 
considered point, uu and a^ are velocity and acceleration of the z-th particle. 

Clearly, Eq. ( TTTT ) should be consistent with Newton's second law governing 
the motion of particles. For a droplet of mass m 8 it reads: 

midLu = m^g + fu. (26) 

Here fu is the force exerted by the air on the z-th droplet. 
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Summing ( ESj ) over n considered particles and dividing the sum by volume 
V we obtain: 

En sr-^n \r~\n o 

y — Plg + *lj Pi - y J b l = y • (/'J 



Force F/ (123) represents the cumulative drag force exerted by the air on 



all droplets in a unit air volume. Comparing Eqs. ( 1271 ) and ( 1TTT) we find that 
Eq. ( TTTT ) can be obtained from Eq. (J2SD only if the following non-linear rela- 
tionship holds: 



En 

v — = pl 



- w + (ur V)u ( 



, (28) 



We will now show that Eq. (J2HD does not generally hold for a flux of particles. 
We will consider the simplest case when all the particles have the same mass, 
rrii = m and pi = nm/V. Eq. (J2HD then reduces to 

ai = ¥ +(ur v)u/ = m- (29) 

For the purposes of our demonstration we divide the particles in the flux 
into two arbitrary types with respect to dynamics: e.g., let us colour these 
types green and red. As the particles do not interact, the flux of green particles 
is not affected by the presence of red particles, and vice versa. Thus, if (J22) 
were generally true, it must apply to the red and green fluxes separately as well 
as to the combined flux. 

Mean velocity u/ and mean acceleration a; of the flow as a whole (red+green) 

are 

ui = aun + Qua, ai = aa a + /3a^ 2 , a + /3 = 1, (30) 

where a and (3 are the relative numbers of red and green particles in the total 
flow, Uik and ( 1251 ) are the mean velocity and acceleration of green (k = 1) 
and red {k = 2) particles, respectively. 

Let relationship flZHD be satisfied separately for both types of particles, i.e. 
the following is true for k = 1, 2: 

a^ + ^.VW,^,, ,3!) 

Putting (15TT) and (13UT) into (12~9T) , we find that (j22D only holds for the combined 
flow if 

([Uii-Ui 2 ] • V)(u n -u Z2 ) = 0. (32) 

This is true if either un = U/ 2 or Vu/i = Vu^ 2 . These conditions apply only 
when all particles possess the same velocity. 

Such conditions do not hold in a precipitating atmosphere where some 
droplets are already falling at near terminal velocity and others have only begun 
to form and have near zero motion with respect to the surrounding air. Droplets 
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of the first type have approximately zero mean acceleration an ~ <C g, while 
droplets of the second type have acceleration = g- Condensation produces 
droplets of the second type only and thus does not change the total momentum 
per unit volume. Moreover, at their instant of forming having zero velocity with 
respect to the air, the newly formed droplets do not impose any force on the 
air. Therefore, F a cannot depend on condensation rate and formula (J2J) must 
be incorrect. 

Finally, let us consider the situation with condensate particles having one 
and the same velocity. Condensation occurs at height z = h only. Consider a 
droplet that arises from condensation at time t = and starts falling. Eq. ( I2B1 ) 
in this case has the form (we omit indexes / and i): 

a = -77 = ~9 + — • (33) 
at m 

Location Z of the droplet on its trajectory depends on time t as 

t <iZ dZ 

u(t')dt\ -i7 = u(t), dt = ——. (34) 
at U \Z J ) 

Changing variables from ttoZ in ( 1551 ) with use of ( 1541 ) we can re- write Eq. ( 1551 ) 

as 

du f d u 2 f 

a = M lz = - 9+ ™' lzY = - g+ ™ (35) 

Eq. ( 1551 follows from the linear equation ( 1551 ) and represents the law of energy 
conservation. Non-linearity of Eq. ( 1551 ) appeared when t was replaced byy z 
in ( 1551 ). Coordinate Z = Z(t) in ( 1551 describes the position of droplet on its 
trajectory and is unambiguously determined by the dynamics of the droplet 
( 1551 ) and time t. 

If force / acting on a droplet at point Z(t) = z is time-independent, all 
droplets forming at z = h will have one and the same velocity u(z) at point 
z. There will be a stationary spatial distribution of particle velocity u(z) (even 
if condensation rate and, hence, particle density vary in time). We can see 
that in this case Eq. (J22D is fulfilled taking into account that du/dt = 0. On 
the other hand, if force f(z) depends explicitly on time, then particles formed 
at different times can have different velocities at point z. In this case, particles 
with different velocities can be present at any point z, which, as we have shown, 
makes Eq. ( 1291 ) invalid. 



7 Conclusions 



Recently ICotton et all [201 1| reflected upon the fundamental equations describ- 
ing moist at mospheric dyna mics. They concluded that the reactive motion term 
descri bed by Bannon |2002 | was of potential importance and warranted further 
study. ICotton et all [201 1| offered various graphical illustrations explaining how 
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the phenomenon is manifested. Here we have shown that this reactive motion 
term was based on a physical error. It does not exist. 
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